Correlation effects in quasi one dimensional electron wires 
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We explore the role of electron correlation in quasi one dimensional quantum wires as the range 
of the interaction potential is changed and their thickness is varied by performing exact quantum 
Monte Carlo simulations at various electronic densities. In the case of unscreened interactions with 
a long range 1/x tail there is a crossover from a liquid to a quasi Wigner crystal state as the density 
decreases. When this interaction is screened, quasi long range order is prevented from forming, 
although a significant correlation with 4&F periodicity is still present at low densities. At even 
lower electron concentration, exchange is suppressed and the electrons behave like spinless fermions. 
Finally, we study the effect of electron correlations in the double quantum wire experiment [Steinberg 
et al., Phys. Rev. B 77, 113307 (2006)], by introducing an accurate model for the screening in the 
experiment and explicitly including the hnite length of the system in our simulations. We find 
that decreasing the electron density drives the system from a liquid to a state with quite strong 
4/cf correlations. This crossover takes place around 22 nm~ , near the density where the electron 
localization occurs in the experiment. The charge and spin velocities are also in good agreement 
with the experimental findings in the proximity of the crossover. We argue that correlation effects 
play an important role at the onset of the localization transition. 

PACS numbers: PACS numbers: 73.21.Hb,71.45.Gm,71.10.Pm 



I. INTRODUCTION 

It is well known that the effect of interactions in 
quasi one dimensional (Q1D) systems of electrons, usu- 
ally called "quantum wires" , is enhanced compared to 
higher dimensional systems. There are universal prop- 
erties described by the Luttinger liquid paradigm, the 
effective low energy theory which applies for strictly 
ID models^^i 3 -^ such as spin-charge separation, charge 
localization, and conductance quantization. However, 
the microscopic details, such as the width and type of 
the transverse confinement or the distance and shape of 
neighboring screening media, can have a large impact on 
the properties of the Q1D systems, as they are very sen- 
sitive to the effective interaction. These systems can be 
realized in semiconductor structures, where there are ele- 
gant experimental studies j&L&i y&lItlS and it is essential 
to describe the system accurately for a realistic compar- 
ison of theory and experiment. 

In this paper we study how the thickness, finite size, 
and screening affect the phase boundaries of some uni- 
versal features with a particular emphasis on the charge 
localization and spin properties. We address the issue 
of how the electron correlation depends upon the mi- 
croscopic details parametrized in the interaction using 
ground state quantum Monte Carlo (QMC) methods 13 
such as diffusion Monte Carlo 14 ' 15 , and its lattice regu- 
larized version^ which are ideal numerical tools to study 
Q1D systems, since they provide exact results in one di- 
mension. Previous QMC studies regarded the determina- 
tion of the LL parameters for a Q1D system with screened 
interactions,^ and the ground state properties of a model 
with a long-range Coulomb potential^ Here we compare 
various model interactions in a unified picture, with the 



final goal of quantifying the role of correlation in the lo- 
calization transition found by Steinberg et al.& 

With this aim it is particularly important to include 
both the effects of the long-range Coulomb potential and 
the consequences of its screening in the interactions. In 
dimensions larger than one, the 1 jx tail of the Coulomb 
pairwise potential leads to a Wigner crystal phase of the 
homogeneous gas at low densities, when the potential 
energy dominates over the kinetic contribution ^ 9 : 20 : 21 In 
this regime the spin exchange drops to an exponentially 
small valued 2 - as the overlap between unlike spin parti- 
cles is exponentially suppressed by the localization of the 
electrons. Therefore, one of the signatures of the liquid- 
to-crystal transition is a decrease in the spin stiffness. 
However, it is possible that other spin and charge phases 
could exist in between. 

In ID the situation is radically different. It is well 
known that a Luttinger liquid (LL) with 1 jx interactions 
exhibits slowly decaying charge-charge correlations, but 
no true long range order, as the quantum fluctuations are 
stronger in lower dimensions. 3 Nonetheless there should 
be a crossover from a high density liquid to a low den- 
sity regime with quasi long-range charge order also called 
a "fluctuating Wigner crystal" ^ However the LL the- 
ory does not predict where the crossover happens, as the 
correlation function parameters are not universal but de- 
pend on the details of the interaction. Also the inter- 
play between charge and spin is quite unclear. Indeed 
the LL parameters depend on the effective one dimen- 
sional potential in a non trivial way. For instance, the 
spin properties are strongly affected by its short-range 
behavior which includes the effect of the thickness, as 
the transverse dimension can effectively tune the spin 
exchange. On the other hand, the quasi order of the 
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charge degrees of freedom is stabilized by the long range 
tail. The relative importance of the short versus long 
range correlations is set by the microscopic model of the 
system. Recently Fogle r 24 i 25 proposed that a correlated 
state with very small spin exchange exists for ultrathin 
wires at densities between the liquid and the quasi Wigner 
crystal phase. In the limit where the short-range part can 
be effectively described by an infinite repulsive contact 
interaction, this state can be related to a noninteracting 
spinlcss Fermi system, as in the Tonks- Girardeau gas>2£ 
In a one dimensional system of fermions, the coexistence 
of strong short-range repulsions and very long-range in- 
teractions leads to a peculiar state, which Fogler termed 
a Coulomb- Tonks gas. 

In previous theoretical work, quantitatively accurate 
studies of the liquid-to-crystal one dimensional crossover 
have been carried out only for inhomogeneous systems, 
with longitudinal extension controlled by an external 
confinement, and where the finite (and very small) size 
allows one to solve the problem by means of exact 
diagonalization . 27 ' 28 ! 29 However, the broken translational 
symmetry leads to quite different properties, particularly 
in the charge and spin density profiles. 

From the experimental side, technological advances in 
the preparation of cleaved edge overgrowth samples have 
enabled tunneling measurements between two high mo- 
bility parallel wires, which probe striking features like 
the spin-charge separations in the excitation spectral In 
a recent extension of the tunneling experiments, Stein- 
berg et al& applied a gate to the upper wire in order 
to tune its electron density by charge depletion. Below 
a critical threshold, measurements revealed a dramatic 
transition which can be interpreted as the onset of lo- 
calization in the wire. Although it is believed that the 
transition is mainly driven by electron-electron interac- 
tion effects as the liquid phase is in a ballistic regime, 
so far there is no agreement between the critical density 
predicted by theory and the actual experimental value. 
Some features of the experiment, like the fringes of the 
differential conductance in the liquid phase and the first 
two peaks of the tunneling current in the localized phase, 
have been explained in an independent particle picture^ 
and at the mean field level^ The LL theory has been 
applied to describe the general features of the tunneling 
current in the Wigner state by assuming a spin incoher- 
ent regime. 31 However it is unclear at which density the 
spin degrees of freedom become incoherent, as the spin 
velocity measured in the experiment is in disagreement 
with previous numerical estimates. Another open issue 
is related to the fact that only a small fraction of elec- 
trons take part in the localized state. It is clear that an 
accurate microscopic description of the experimental sit- 
uation is necessary to account for all these features. In 
our study we include the most important details such as 
an accurate screening and the effect of the finite size of 
the wire to correctly describe and understand the physics 
underlying the experiment. 

Throughout the paper we use units of the effective 



Bohr radius Oq = " \ for length and the effective Ry- 
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dberg Ryd* = ^f— r for energy where e is the dielectric 
constant of the embedding medium and m* is the effec- 
tive electron mass. 

The paper is organized as follows: in Sec. [H] we present 
results for a quasi one dimensional electron gas (1DEG) 
with long range (1/x) interactions for different densities 
p = ^- and thicknesses. We carefully study the liquid- 
to-quasi-crystal crossover by varying the Wigner-Seitz ra- 
dius r s which sets the relative importance of the kinetic 
energy and the interaction. By means of QMC tech- 
niques, we also study the charge compressibility and the 
spin susceptibility in order to analyze the interplay be- 
tween the charge and spin properties of the wire. The 
effect of the wire's thickness on the crossover and the 
spin properties is taken into account by performing sim- 
ulations with three different wire widths. In Sec. IIIII we 
compare the unscreened 1/x potential with an interac- 
tion screened by a metallic plane. In Sec. IIVI we inter- 
pret the localization transition found in the series of two 
wire tunneling experiments&^&^ i 10 ' 1 ^ 12 by studying the 
evolution of liquid-to-crystal correlations in a finite wire 
with interactions effectively screened by another parallel 
wire. We make a comparison between the finite system 
and the corresponding homogeneous infinite system in- 
teracting with the same potential. We also show the 
agreement between our model and the experiment. Fi- 
nally in Sec. [V] we summarize our results and comment 
on possible refinements to our calculations. 



II. UNSCREENED COULOMB INTERACTIONS 

We study a system of electrons interacting via the 
Coulomb (1/x) potential which are confined to one di- 
mension by a harmonic potential in the transverse direc- 
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tion V(r±) = ^y, where b tunes the thickness of the 
wire. This system was previously studied using QMC by 
Casula et al. and here we follow the conventions used in 
that workji& We integrate over the transverse degrees of 
freedom, which is a good approximation when the density 
of electrons in the wire is low (r s 3> nb/4), and hence the 
longitudinal energy scale is small compared to the excita- 
tion energies related to the perpendicular motion. This 
integration yields an effective one dimensional interac- 
tion: Vb(x) = ^ cxp (jg^j erfc > which has a long 

range 1/x tail. The thickness b of the wire controls the 
short-range behavior of the potential, which is finite at 
the origin (V(0) — ^pn/b). Since the crossover between 
the short and long range behavior is at x ~ b, for smaller 
b the repulsion is stronger as the particles approach each 
other. 

In this work we have chosen to study three different 
thicknesses, b = 1, 0.1, and 0.0001. The first two values 
correspond to typical experimental thicknesses for semi- 
conductor quantum wires, whereas the last one is chosen 
to explore the ultrathin limit as studied analytically by 
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Fogle r 24 i 25 and experimentally realized in carbon nan- 
otubes placed on SrTiC>3 substratesi 32 ' 33 For each value 

of b, the density in the wire ( 2r~J ^ s var i e d, allowing for 
the interaction strength to change and the corresponding 
ground state properties are computed. 

We use diffusion Monte Carlo (DMC) and lattice reg- 
ularized diffusion Monte Carlo (LRDMC) 1 ^ methods to 
project the initial variational ansatz to the ground state 
|^o)- These methods are particularly suited to the sim- 
ulation of one dimensional fermions. Indeed, the well 
known "sign problem" does not affect these calculations 
as the nodes are fully determined by the points of coin- 
cidence between the electrons and therefore are exactly 
included in the trial wave function \^t}- Since the fi- 
nal DMC or LRDMC distribution is the product of the 
true ground state and the trial wave function, some ob- 
servables such as the density and the structure factor are 
determined using the forward walking techniqu o 34 ' 35 in 
order to generate unbiased expectation values. 

We simulate an unpolarized wire with electrons sub- 
ject to periodic boundary conditions (PBC). The trial 
wave function is written in the Slater- Jastrow form 



(1) 



where the Slater determinants for up and down spin elec- 
trons read 



l<i<j<N" 



(2) 



with G = 2ir/L, and L = 2r„N the length of the simu- 
lation cell. We follow Ref. [30 to determine the Jastrow 
function u(x). Its Fourier components are 



2pu(k) = -Soik)-' + V^o(fc)- 2 + 2pV b (k)/k 2 , (3) 

withS (fc) = (k/2k F )d(2k F -k)+9(k-2k F ) the structure 
factor of a noninteracting 1DEG, p — j^- the density, and 

Vb(k) the Fourier transform of Vb(x). To reduce the finite 
size effects in our simulation we use the Ewald technique 
to sum our potential as discussed in detail in Ref. LL8]. 
This approach has been used to study the infinite wire 
with the long range potential Vb(x), and also the screened 
potentials described in the next Sections. In the latter 
case, the sum over the images has been done numerically 
in the real space as the potentials have a shorter range. 

To reveal the presence of charge ordering in the sys- 
tem, we first analyze the static structure factor S(k) = 

(p(—k)p(k)) , where p(k) = Ylj e are the Fourier 
components of the electron density. At high density the 
structure factor is very similar to the mean spherical ap- 
proximation (MSA)2I prediction SmsaW = S (k)/(1 + 
2pu(k)So(k)) as expected (see Fig. [IJ, since in the limit 
r s — > the MSA becomes exact i 1 ^ Specifically, there is 
no peak at 4kp up to r s — 0.5 (r s = 0.2) for b = 0.1 



(b = 0.0001), namely there are no correlations with the 
mean interparticle spacing (Fig. [2]). As the density de- 
creases, a peak develops at Ak F . This peak is a necessary 
feature for a one dimensional quasi Wigner crystal and it 
is absent in the MSA prediction which has no structure 
at ikp. For b — 0.1 we carried out simulations with up 
to 450 particles for r s — 0.5 and r s — 0.75, to check the 
convergence of the S(k) in the liquid regime close to the 
onset of the Ak F charge correlations (Fig. [2]). 

The scaling of the height of the Ak F peak of S(k) 
with the number of particles (reported in Fig. [3] for 
b = 0.1) highlights the features of a liquid-to-quasi- 
crystal crossover. When the peak is absent there is no 
significant dependence of the S(<ikp) value as a func- 
tion of system size, however where there is a peak in the 
structure factor at ikp, its scaling is sub- linear, signal- 
ing a quasi-long range order (linear scaling would indi- 
cate a true Wigner crystal). The points in Fig. [3] are fit 
very well by a functional form obtained from the charge- 
charge correlation function 3 ^ derived by Schula^ 3 - in the 
LL framework with long range interactions, 



da; exp(— iAkpx) (p(0)p(x)) = aL exp(— 4c-\/log L)+b, 

where we explicitly include the dependence on the sys- 
tem size L by taking the Fourier transform over the sim- 
ulation cell. The short-distance cutoff cq is introduced 
because the LL theory provides only the asymptotic be- 
havior for (p(Q)p(x)) . Further logarithmic corrections 
could be included 3 ^ in Eq.[4j but we take just the leading 
order expansion, which should be the most relevant for 
the system sizes computed here. One would need much 
larger systems which are beyond our current numerical 
capabilities to resolve further corrections. The bosoniza- 
tion formalism gives a parameter dependent scaling for 
the Ak F component of (p(0)p(x)) which is left undeter- 
mined in the LL theory, and depends on the details of 
the interaction. At high densities there is no peak in the 
structure factor and the electron gas is liquid (Fig. 
Consequently, there is no finite size dependence at 4kp 
and the parameter a undetermined in the LL theory is 
zero. 

We also determine the charge compressibility \ p and 
the spin susceptibility \a of the electron gas by using 
two techniques in our calculations. The first is to ap- 
ply the definition of those quantities as the reciprocal of 
the second derivative of the total energy with respect to 
the density r s or the polarization Q. Our QMC calcula- 
tions provide measurements of the total energy, so these 
derivatives can be taken by fitting our data with a suit- 
able functional form. The error in such a determination 
comes from both the statistical uncertainty in the calcu- 
lations and the constraint represented by the choice of the 
fitting function. As a technical detail, it is also necessary 
to extrapolate the energy to the thermodynamic limit 
which can be a costly proposition. Moreover a functional 
form which includes the dependence on both density and 
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FIG. 1: (Color online) Structure factor for b = 0.1 (upper 
panel), and b = 0.0001 (lower panel), computed for a system 
with 78 electrons. The QMC (points) and MSA (solid lines) 
structure factors are reported for different densities (r 3 ). Also 
the noninteracting spinless fermion (NSF) structure factor is 
drawn (solid black line) for comparison. 



polarization has not been provided yet4£ Nevertheless, 
we use the parametrization in Ref. Il8l which holds for a 
system with fixed polarization £ and depends only on r s 
to compute the charge compressibility and validate the 
second method to evaluate \ p and \<j- 

The other method we use to compute these quantities 
is to calculate the momentum resolved excitation energies 
of the system, and exploit the sum rules which exactly re- 
late the collective modes of the long wavelength spectrum 
with Xp and x<j- Gold and Calmels^i found that 



uj p (k^0)=v F \k\J PF V{k^0) + ^, (5) 

V *p 



uAk — ► 0) = v F \k\. — , 

V x<y 



(6) 



where cu p (k) (oj a (k)) is the energy of the lowest charge 
(spin) excitation with momentum k, p F is the density of 
states of the free electron gas at the Fermi energy, and 
Xo = 16rj?/7r 2 is its compressibility. 




FIG. 2: (Color online) Detail for the structure factor near 
Ak F for b = 0.1, computed for N = 182 and N = 450 at two 
densities (r, = 0.5 and r s = 0.75) in the proximity of the 
crossover from a liquid to a quasi-crystal. 
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FIG. 3: (Color online) Scaling of the 4k F component of the 
structure factor with respect to the number of particles. The 
scaling is reported for various densities with b — 0.1. The 
lines are the best fit of the function in Eq. U given by the LL 
theory. 



In order to find out the lowest energy states of a given 
momentum k we employ a method proposed by Ceperley 
and Bemuds, which is a generalization of the transient 
estimate used in the projection Monte Carlo (DMC or 
LRDMC) framework. This method is based on the idea 
that it is possible to compute the excitation spectrum 
of a system in a direct and variational way by project- 
ing the initial basis functions to their lowest energy state 
with the given symmetry. In our case the basis set is the 
Feynman ansatz-S i.e. p(k)\^o) Vfc for the charge exci- 
tations and cr(fc)|\f r o) Vfc for the spin excitations, where 
a{k) = J2j J2a crelkTj is the Fourier transform of the 
spin density. In the following we assume to work with 
the charge excitations, but the same applies for a(k). 
Since the basis set is orthogonal, the method in Ref. |42| 
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is greatly simplified, as every k component is decoupled. 
For each k, we have to calculate 

(M> |p(fc,r)gp(-fc ; o)|^ ) _ Ej44^!!^ m 

{9 \p{k,T)p(-k,0)\* ) Ei4e-(4-^o) ' [) 

where p{k,r) is written in the Heisenberg representa- 
tion with imaginary time evolution, VP^) is the i'th ex- 
cited state with momentum k, el is its energy, A\ — 

l 1 2 

I (^fe|p(~^)l*o) | is the spectral weight of the eigenvalue 
expansion, and Eq is the ground state energy. For large 
r the ratio in the above Equation will converge to the 
lowest energy e° of a given k, provided A® is non zero. 
Another limitation is given by the exponentially small 
denominator, which will exponentially increase the sta- 
tistical noise of the estimate as the projection time in- 
creases. Both the numerator and denominator in Eq. [7| 
are evaluated by means of the forward walkin g 34 ' 35 pro- 
cedure based on the DMC or LRDMC sampling. Indeed, 
for large enough r the LHS of Eq. [7] can be rewritten as 

J dndr 2 p(-k)G( ri ,r 2 ,T)E L (k, r 2 )p(k)P(r 2 ) 

/dndra p(-k)G(r u r 2 ,r)p(k)P(r 2 ) ' U 

where E L (k,r) = "jffiyjfy^ is the local energy of 
p(k)\^ T ), P(r) = ^ T {r)^o(r) is the QMC mixed dis- 
tribution, and G(ri,r 2 ,r) = ^ T {ri){r 1 \e-' rH |r 2 >/*r(r 2 ) 
is the importance sampled Green's function. 

Because the excitation energies uo{k) — e" — Eq are 
computed relative to the ground state energy Eq, there 
is a cancellation of errors since the sample generated to 
compute Eq and e° is the same. Therefore a modest size 
calculation is enough to get converged energies. The con- 
vergence with the propagation time can be more difficult 
to obtain. However, for the long wavelengths p(k) \$q) 
is a good approximation to the lowest excited state with 
momentum k and the energies can be determined easily 
with a short projection time r. When the small k range 
of energies is fit to the form in Eqs. [5] and [SJ Xp an d 
Xa- are determined. The results for the charge compress- 
ibility obtained with this method agree with the second 
derivatives of the total energy in all cases we have made 
the comparison, as is shown in Fig. 0J 

The knowledge of Xp an d X& can shed more light on 
the properties of the liquid-to-quasi-crystal crossover. By 
looking at the charge compressibility (Fig. 0|, it is ap- 
parent that the role of the electron correlation is be- 
coming increasingly important in the proximity of the 
crossover, where there is significant discrepancy between 
the Hartree-Fock (HF) and QMC values of Xp- I n P ar ~ 
ticular, the correlation makes the system softer than the 
HF, which is consistent with a more pronounced localiza- 
tion of the electrons. At even lower densities the charge 
compressibility of the unpolarized system is approach- 
ing that of a fully polarized (or spinless fermion) gas. 
The difference between the two is going exponentially to 
zero, and they almost overlap for r s > 4 (with b = 0.1). 
This means that the energy of the spin excitations is get- 
ting smaller and smaller as the density decreases. This 



feature is revealed by the inverse spin susceptibility Xa- 
The Xo/x<t ratio is plotted in Fig. [5] This value be- 
comes exponentially small at low densities, where it is 
difficult to get a statistically accurate QMC estimate, 
since the sampling of the spin is "frozen" by the pres- 
ence of quasi nodes (pseudo nodes) between unlike spin 
electrons^ The strong interaction makes the electrons 
to repel each other at short-range, and the corresponding 
wave function is very small at the coalescence points of 
electrons with opposite spin. Consequently the spin flip 
rate in the QMC sampling becomes small, and the effi- 
ciency decreases. However the charge properties do not 
seem to be affected by this slowing-down. The physical 
reason for the quasi nodes will become even more appar- 
ent in Sec. IIIIl when we will discuss the Tonks- Girardeau 
physics of the screened wire. 

In the low density regime where exact Monte Carlo 
sampling becomes difficult the WKB approximation is 
useful for determining the dynamical properties of the 
electron gas. Following the example of Matveev^ 4 - we use 
the WKB approximation to determine the rate at which 
two electrons exchange by calculating the energy barrier 
that they must overcome. Although fluctuations prevent 
the formation of a Wigner crystal, the equilibrium posi- 
tions of the electrons are assumed to be equally spaced 
with periodicity 2r s . Central to the accuracy of this ap- 
proximation is the fact that at low densities the tunnel- 
ing is dominated by the effect of the potential, and the 
statistics can be ignored. Furthermore, all electrons are 
treated as uncorrelated except for a single pair which is 
allowed to exchange. In contrast to Matveev's approach 
we assume that the other electrons are distributed about 
their equilibrium positions according to the harmonic 
approximation with a Gaussian spread instead of being 
fixed delta function point particles. Taking the initial 
positions of the two exchanging electrons to be at x = 
and x = 2r s , they feel a static potential given by 

/oo 
p(y)V(x - 2nr s + y)dy, (9) 

where p(y) = \fapn exp(— ay 2 ) is the equilibrium charge 
density of the non exchanging electrons and V(x) is 
the interparticle potential. The harmonic approximation 

gives a = J m^-^^-, where W(x) is the potential at 
a given lattice site due to an infinite array of electrons 
spaced as 2r s . 

At low densities the electrons behave as a spin chain 
obeying the Heisenberg Hamiltonian were the spin flips 
are mediated by an exchange of nearest neighbor elec- 
trons, so the spin susceptibility can be determined from 
the energy barrier computed within the WKB approxi- 
mation by analogy with the Heisenberg Hamiltonian in 
ID as shown by Matveevi 4 ^ The spin velocity of the 
equivalent Heisenberg spin chain can be found from the 
Bethe ansatz solutio n 45 ' 46 , yielding v a = irJr s /2 where 
J is the size of the energy barrier in the WKB approxi- 
mation. This gives the susceptibility through Eq. [5] 
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Where the density is large enough that QMC reliably 
samples the spin exchanges the spin susceptibility com- 
puted using the forward walking techniques agrees well 
with the WKB estimate only after the smearing of the 
electron sites given by the harmonic approximation. It is 
therefore important to use the potential in Eq. [9] to have 
an accurate estimate of the exchange at intermediate den- 
sities. This agreement and the fact that the dynamical 
many-body corrections to the WKB estimate are very 
small at low density 7 justify the use of WKB for dilute 
systems where it is difficult to extract information from 
the QMC calculations. In addition, the exponential de- 
cay of v a versus -y/rj obtained in this way is in agreement 
with previous result o 44 ' 47 ' 48 for potentials where they can 
be compared. 

Fig. [5] summarizes our findings for the unscreened wire. 
The liquid-to-quasi-crystal crossover is shifted to higher 
densities for thinner wires, while the spin susceptibility 
is always significantly different from zero in the crossover 
region for the values of the confinement taken into ac- 
count. The smallest b we studied (b = 0.0001) cor- 
responds to one of the thinnest confinements realized 
experimentally^^ The spin exchange is still sizable in 
the crossover region due to the not-so-long localization 
length of the electrons and not-so-thin width of the wire. 
Therefore, in our study we did not find any signature 
of the Coulomb- Tonks gas phase in between the liquid 
and quasi Wigner crystal, which was claimed by Fogler 
for ultrathin wires. 25 However, the structure factor plot- 
ted in Fig. Q] reveals the tendency for electrons to ap- 
proach the noninteracting spinless fermion behavior (the 
limit where the Coulomb- Tonks gas picture holds) as the 
wire width decreases. The fundamental difference with 
respect to the noninteracting spinless picture is the pro- 
nounced peak at Akp, which characterizes the Coulomb 
long-range interactions at low density. 




FIG. 4: (Color online) Inverse charge compressibility Xo/Xp 
of the unpolarized and fully polarized wire for b = 0.1, with 
both screened (R — 1) and unscreened interactions. Also 
the HF (dashed back line) charge compressibility is reported 
for the unpolarized wire. The solid lines are obtained from 
the second derivative of the energy parametrization, while 
the points are evaluated through the charge excitations as 
explained in Sec. [II] 
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III. SCREENED INTERACTIONS IN GATED 
WIRES 



The primary interest of this paper is to model a quan- 
tum wire formed in a semiconducting nanodevices. In 
that case there is almost always a metallic gate that 
screens the long range (1/x) potential. To see the changes 
that such a gate would cause, we introduce a perfectly 
conducting metal plane parallel to the wire located a dis- 
tance R away. Using the electrostatic method of images 
the potential is constructed by assuming that a wire is 
placed at a distance 2R from the original one with the 
same particle distribution but opposite sign. The equa- 



FIG. 5: (Color online) Inverse spin susceptibility xo/x<? f° r 
different thicknesses. The dependence on r s is shown. The 
points are the QMC calculations, while the lines are the WKB 
estimates. The arrows indicate the liquid-to-quasi-Wigner- 
crystal crossover. 



tion for this potential is 



V(x) 



drdr 



drdf' 



, Pb{r)p b {r') 
y/ (r — r') 2 + x 2 
Pb(f)Pb(r') 



\r-r' -2R) 2 + x 2 
= V b (x)-V mt (x,R) (10) 

where r and f' are transverse vectors, Pb(r) — 
b\/2n ex P ( ~itp) * s ^ c ground state charge distribu- 



7 



tion of a two dimensional harmonic oscillator with the 
wire's confining potential: V w i IC (r) = The first in- 
tegral gives the effective unscreened inter-particle poten- 
tial Vb(x) described in the previous section and the sec- 
ond one is the potential due to the image charge on the 
screening wire: Vi n t(x, R). 

The quasi Wigner crystal correlations derived by 
Schulz 23 apply only when the interaction is long range 
{!/%)■ In the case of the screened interaction above the 
potential decays as 4R 2 /x 3 at large distances, so a simple 
scaling argument shows that the Wigner crystal correla- 
tions should be absent at very low densities. Indeed, if 
r s > 8R 2 /it the typical kinetic energy of the electrons, 
the Fermi energy Ef , is larger than the potential energy 
computed at the mean interparticle distance (2r s ). At 
these low densities Matveev^ 4 - has pointed out that it is 
possible to map the screened short-range interaction into 
a repulsive contact potential 

V{x) = US(x), (11) 

where the constant U is chosen so the delta function 
potential and the screened one have equal transmission 
coefficients. On the other hand, in the density range 
1 <C r s < 8R 2 /ir the 1/x shoulder of the potential can 
induce Akp correlations, which are strong but not strong 
enough to stabilize any sort of quasi-order. Calculations 
of the finite size scaling of the 4/cf peak of the structure 
factor for b = 0.1 and R — 200 show the saturation of its 
height for N > 100, and so demonstrate the absence of 
the quasi Wigner crystal correlations when screening is 
introduced despite quite a large distance to the metallic 
gate (FigHJJ). Only in the limit of R — > oo does one re- 
cover the unscreened potential and the possibility for a 
quasi long-range charge order. 
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FIG. 6: (Color online) Scaling of the 4&f component of the 
structure factor with respect to the number of particles for 
b = 0.1, and r 3 — 4. For comparison, the scaling is reported 
for the unscreened V\> interaction, and the screened potential 
in Eg. niTl with R = 200. 

The lack of the quasi Wigner crystal state does not 
change the crossover to the spinless fermion physics 



present in the unscreened system. Even though the quasi 
long-range charge order is absent at low densities, the 
charge compressibility still approaches that of a gas of 
spinless fermions as the density decreases. This approach 
can seen in Fig. [4j It is therefore clear that the spin 
crossover does not depend on the long-range correlations. 
In fact, this crossover can be reproduced by a system of 
electrons interacting via the delta function interaction in 
Eq . ITTI where the constant U is large, an interaction that 
has no long-range piece whatsoever. 

The low density limit with screened interactions is par- 
ticularly interesting as the screening introduces a new 
feature. At low densities the electron-electron repulsion 
at short range makes exchanges between electrons virtu- 
ally impossible, corresponding to the limit U —too. As a 
result for the ultrathin wire with strong screening (b <C 1 
and r s 8R 2 /ir), the mapping of the interaction to the 
potential in Eq. 12 becomes exact. In this situation not 
only do the electrons behave as spinless fermions, but the 
charge velocity approaches that of noninteracting spin- 
less fermions (v p = 2vf)- This is analogous to the case 
of bosons with infinite repulsive contact interactions, (or 
impenetrable particles) where the system can be mapped 
into a noninteracting Fermi gasi2£ The impenetrable Bose 
system is often called a Tonks- Girardeau gas. In our case 
the situation is analogous, namely the fermions become 
impenetrable due to an effective infinite contact repul- 
sion, and so they behave as they were noninteracting and 
spinless. We refer to this behavior as Tonks-Girardeau 
regime. One of its features is the presence of nodes in 
the wave function at the coalescence of unlike spin pairs. 
This is the extreme case when the pseudo nodes that 
complicate the ergodicity of Monte Carlo calculations at 
low density as reported in Sec. HT1 become actual nodes. 

While this effect has been discussed in the 
literature ; 24 ' 25 ! 44 our work provides quantitative predic- 
tions for the onset of the noninteracting spinless behavior. 
Fig. [7] shows the charge velocity in the limit of low den- 
sity for different values of the screening in the thinnest 
wire we studied (b = 0.0001). We found that in order for 
the Tonks-Girardeau behavior to manifest itself, the dis- 
tance to the gate R must be less than 0.1 and the density 
must be lower than r s = 1. For R larger than 0.1, at low 
density the charge velocity does not converge to the non- 
interacting spinless fermion limit (2vf), but saturates at 
a larger value. 

It is possible to see the transition of the screened elec- 
tron gas to the noninteracting spinless fermion behavior 
more directly by analyzing the static structure factor, as 
was done in the unscreened case. In Fig. [5J the S(k) is 
plotted at different densities for the ultrathin wire with 
b = 0.0001 and gate located at R = 0.1 from the wire. 
Contrary to the case of the unscreened wire (Fig. [1] lower 
panel), at low densities the peak at 4fcp is absent and 
the structure factor approaches that for noninteracting 
spinless fermions quite closely. Notice that at the same 
time the charge velocity approaches the value of 2vf (see 
Fig. ED. 
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FIG. 7: (Color online) Asymptotic large r s values of the 
charge velocity in units of vf vs. inverse screening length 
for ultrathin wire (b = 0.0001) from R = 0.05 to R = 5. In 
the inset we report the full dependence of the charge velocities 
on r s at different R. 
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FIG. 8: (Color online) Static structure factor for the screened 
wire with b — 0.0001 and R — 0.1, plotted for three values of 
the density, r s = 0.2,0.4 and 0.8. The solid lines correspond to 
the MSA prediction for each density, and the black line is the 
structure factor for noninteracting spinless fermions (NSF). 



The same study was repeated for the wire with b = 
0.1. Here the short-range behavior of the potential is 
much less repulsive than in the b = 0.0001 case and the 
same value R for the screening. The result of this is that 
the charge velocity does not converge to 2vf even for a 
gate as close as R — 0.1, which equals the width of the 
wire and thus represents the geometric limit of validity 
for the uncorrelated inter-wire interaction. Therefore for 
b = 0.1 and thicker wires, whose widths are realizable 
in semiconducting nanostructures, we did not find the 
Tonks-Girardeau behavior in our calculations. 



IV. LOCALIZATION TRANSITION IN TWO 
PARALLEL GaAs WIRES 

Quasi one dimensional systems can be realized in 
GaAs/AlGaAs heterostructures by means of various 
techniques. One such technique being cleaved edge over- 
growth, which has been applied recently to build an ex- 
perimental setup with two parallel wires so that it is pos- 
sible to observe momentum resolved tunneling from one 
to the other l £iL5i2 a 10 i 11 i 12 In this series of experiments 
both the energy of the tunneling electrons and their mo- 
mentum could be tuned by changing the relative chemical 
potential and the applied magnetic field. This setup al- 
lows the dispersion relations of each wire to be probed 
in a quite straightforward manner. Steinberg et al£ fur- 
ther explored how this tunneling is affected by a gate 
that depletes the density of the electrons in the upper 
wire. They found that as the density is decreased there 
is a marked transition in the tunneling interpreted as a 
transition from a liquid to a localized state. 

In the experiment, the center to center distance be- 
tween the two wires is R = 31nm. The upper wire is 
2 \xm long and 20 nm wide. It is the probe to study the 
electron localization. The electrons tunnel from the lower 
wire, which has a width of 30 nm and is taken to be in- 
finitely long. This is also a screening medium for the up- 
per wire. The system is fabricated out of GaAs for which 
e = 13.1 and the effective electron mass is m* = 0.067 m e . 
This gives an effective Bohr radius aS = £ " ■> w 10 nm. 
For the experiment in question the electron density in 
the lower wire is around 60/im _1 , which corresponds to 
r s = 0.83 in units, while in the upper wire the density 
is varied by tuning the gate voltage Vg ■ The effect of Vg 
on the lower wire is very small^ and can be neglected. 

The results presented in the previous section offer an 
avenue to explore the role of the electron correlation in 
the transition observed in the experiment. As the den- 
sity in the wire decreases the strength of the potential 
increases relative to the kinetic energy. One effect of this 
increased relative strength is that exchanges between the 
electrons are suppressed, causing the system to crystal- 
lize. To better quantify the importance of this effect in 
the experimental system, in this section we take into ac- 
count a more realistic potential, assuming the electrons 
are screened by the lower wire instead of an infinite 
metallic gate. To construct this interaction we neglect 
the correlation between the wires and treat the screen- 
ing effects coming from the electrons in the lower wire 
within the linear response theory. We write the potential 
in Fourier space 

V{k, R) = V b (k) + Vim [k, R) X (k)V int (k, R), (12) 

where Vb(k) and Vi n t(k, R) defined in Eq.[]I)]are the intra- 
and inter- wire potentials respectively. Vi n t(k, R) is eval- 
uated by assuming that the thickness of the two wires is 
the same (and equal to the upper wire). This significantly 
simplifies the form and the calculation of the inter-wire 
interaction. x(fc) is the static density-density response 
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function of the lower wire, taken in the random phase 
approximation (RPA): 



XRPA{k) 



l-V b ,(k) X o(k) 



(13) 



fc— 2fej 



k+2k F 



is the static response func- 



where xo{k) = ^ In 
tion for a one dimensional noninteracting Fermi gas, and 
b' is the width of the lower wire. The experimental ge- 
ometry sets the parameters in our quasi one dimensional 
interaction V(k,R). The confinement potential for the 
upper wire is chosen so that the electrons are constrained 
to be inside the lOnm thick wire. Specifically, we require 
the radial root mean squared displacement is equal to the 
lithographic thickness yielding b = 0.707(« l/\/2) for the 
upper wire. The choice of confinement also agrees well 
with the experimental observation that a second mode 
becomes populated at n = 80/xra -1 J£ Similarly, the 
lower wire's thickness is given by b' — 1.061(« 1.5/y2)- 
The distance between the wires is R — 3.0, while the 
Fermi momentum in the RPA response function for the 
lower wire is set by the density r s — 0.83. 

Our screened potential in Eq.[l2]is similar to that used 
by Fiete et al^L who chose a perfect metal response func- 
tion which is valid when the screening wire is at very high 
densities. Here we use the RPA which depends on the ex- 
perimental density of the lower wire through the value of 
the Fermi momentum kp. We notice that our screened 
potential equals that in Ref. [HI at k = 2kp and in the 
limit of small k, namely the long-range tail is the same, 
decaying approximately as l/x*. 

We first analyze the homogeneous system and then ex- 
plicitly include a longitudinal confinement in our simula- 
tions to quantify the finite-length impact on the proper- 
ties of the system, and more closely reproduce the experi- 
mental situation. In the homogeneous system of electrons 
interacting via the potential in Eq. [121 we observe the ap- 
pearance of a Akp peak in the S(k) around r s = 2.2. As 
shown in Fig. [51 it is clearly visible for r s > 2.6, whereas 
no peak is discernable for r s 1.9. This crossover is 
similar to that found for long range l/x interactions. 
However, the important difference here is that the quasi 
long-range order is not present in this case. Indeed, we 
have made a systematic study of the scaling with size, 
and the height of the peak converges to a finite value in 
the thermodinamic limit for all densities taken into ac- 
count. This behavior is consistent with the decay of the 
screened interaction which is faster than l/x. 23 There- 
fore, the crossover is between a high-density liquid to one 
with strong Akp correlations, whose onset can be related 
to the transition occurring in the experimental system. 

The above treatment of the upper wire as infinite and 
homogeneous can be improved to resemble the experi- 
ments more closely. In the study of the 1DEG there 
are strong effects due to any perturbation that breaks 
the translational invariance of the system. For instance, 
Tserkovnyak et al. showed that the asymmetry in the 
oscillations of the conductance as a function of the mo- 
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FIG. 9: (Color online) Static structure factor for a homo- 
geneous wire with b = 0.707 interacting with the effective 
potential in Eq. 1121 which includes the screening by another 
homogeneous wire with r s = 0.83, b = 1.061, and R — 3. The 
structure factor is plotted for several values of the upper wire 
density, with r s ranging from 1.7 to 9.4. The calculations 
have been converged to the thermodynamic limit, requiring 
TV = 62 for r s < 3.0 and N — 78 subject to periodic boundary 
conditions for r s = 4.6 and 9.4. 



mentum transferred between the two wires can be ex- 
plained at the WKB level by having a soft confinement 
potential for the upper wire4 In a later paper they accu- 
rately determined the functional form of the longitudinal 
confinement by fitting its parameters to reproduce the 
period of those oscillations as a function of the magnetic 
field applied to the sampled The potential that provided 
a good fit to their data reads 



V(x) 



Ep { 2 4 



(14) 



where Ep is the Fermi energy of the upper wire, and L 
is approximately 1.5 times the lithographic length of the 
upper wire, namely L = 300 in units. 

We used the above potential together with the inter- 
particle potential in Eq. [12]to study the effect of the con- 
finement on the transition. Although in principle diffu- 
sion Monte Carlo yields an unbiased ground state energy 
in one dimension even for a confined system, (the nodes 
being exactly determined by the coalescence conditions 
just as in the infinite homogeneous wire) in practice it is 
necessary to improve the guidance wave function to re- 
duce the variance of our estimates. The Jastrow factor 
used in the homogeneous system (Eq. [3]) is replaced by a 
more sophisticated factor including one-, two-, and three- 
body terms, fully optimized by means of the stochastic 
reconfiguration (SR) algorithm , 49 ! 50 while the Slater part 
is kept the same as in Eq. [U The one-body Jastrow 
exp(Ji) is needed to localize the electrons in the finite 
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system. It reads 



* = E(- 



(3x 



1 1 > 



(15) 



where a is a free parameter and (3 — \/ E u (2 / L) 4 " / b is 
fixed to cancel the contribution of the potential to the 
local energy at the leading order in the large distance ex- 
pansion. The two-body exp(J2) and three-body exp(J3) 
Jastrow factors are given by 



■h 



(Xjj), 



and 



(ia),(jo-'),(ko-") 



U%" {Xij) U 



a a 
3 



(16) 



where the interparticle distance. Since the finite 

system with screened interactions is dominated by short- 
range correlations, we chose u n (x) to have a simple Gaus- 
sian form 



\x) = Sf exp (-£77,7') 



(18) 



with <5^ CT and variational parameters. Energy min- 
imization improves the quality of the variational wave 
function and stabilizes the forward walking estimate^ of 
the expectation values on the DMC projected state. 

Again the static structure factor is determined for dif- 
ferent densities of electrons in the upper wire. In contrast 
to the calculations for the homogeneous system, the den- 
sity of the electrons is not a direct input to the calcu- 
lation. Instead, we control the number of electrons in 
the wire which are then free to relax according to the 
external potential. An average density can be deter- 
mined by considering the locations of the 2k p and Akp 
peaks of the structure factor and comparing their value 
to those of an infinite array of electrons, 2kp = ^fr and 

4fcf = f-, f s being the effective density in the system. 
Using these conventions, the structure factor for several 
different numbers of electrons is plotted in Fig. [10] 

In addition to the formation of a broad peak in the S(k) 
at 4fcp around N — 80, which corresponds to r s = 2.3, 
the density profile n(x) — {^2 i 5(x — Xi)) of the electrons 
also shows a clear cut sign of the transition. At low den- 
sities, the electrons are distributed in order to minimize 
the interparticle repulsion. This leads to N oscillations in 
the density profile of the wire, a configuration also called 
"Wigner molecule" j£! which corresponds to the 4kp peak 
in the S(k). When the density is increased, the number 
of peaks in the density profile is reduced by a factor of 
two, the Pauli exclusion principle between like spin parti- 
cles being the only factor that prevents the electrons form 
crossing each other. At the same time the Akp peak in 
the S(k) disappears and only a 2kp singularity is present. 
The density is plotted in Fig. [11] for half of the wire as 




FIG. 10: (Color online) Static structure factor for a wire as 
in Fig. [9] but with finite length S(k) is plotted for 20, 40, 
60, 70, 80, 90 and 100 electrons. The corresponding effective 
densities f s are reported in the legend. 




FIG. 11: (Color online) Density profile for electrons in the 
finite wire as in Fig. 1101 plotted for half of the wire length. 
N = 20, 40, 60, 70, 80, 90 and 100 are considered. 



the system is symmetric under inversion around its cen- 
ter. This plot also suggests a transition near N = 80. 

Surprisingly, the calculations with the confinement po- 
tential and the infinite wire give very similar structure 
factors in the vicinity of the transition, suggesting that 
the interparticle correlations are not strongly affected by 
the external confinement at those densities. At lower 
densities the peak at Akp is much larger for the homoge- 
neous system because of the limited number of particles 
in the finite wire. Both the infinite and finite wires show a 
transition from a system with 2kp correlations to a state 
where correlations have a 2r s periodicity. The crossover 
occurs around r s — 2.3, which corresponds to the density 
of 22/im -1 in a GaAs heterostructure. This is very close 
to the density found by Steinberg et al. (20/im _1 ) for 
the localization transition in wires where one subband 
is occupied. However it seems that in the experiment 
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the localization involves only few particles (up to 12 in 
the highest density localized state), i.e. only a section of 
the wire takes part in the transition. This is an impor- 
tant difference with respect to our calculations where the 
transition takes place throughout the system in a quite 
homogeneous way. A non homogeneous behavior is found 
at the edge of the wire where the confining potential in 
Eq. [H] turns upward. There the transition happens at 
higher densities, as one can see in Fig. [TTJ This can be 
understood in terms of a local mean field description. 
At the edge of the wire the effective chemical potential 
fj,0 — V(x) is smaller, corresponding locally to a fluid at 
much lower density. 

Apart from these features, we did not find any Wigner 
correlated patch embedded in a liquid-like system, which 
seems to be the experimental outcome. Therefore a more 
detailed analysis of the experimental setup is required 
to understand better the experiment. For instance, one 
of the top metallic gates used to tune the upper wire 
density could induce a plateau in the external potential, 
nucleating a Wigner region as suggested by Mueller.— On 
the other hand, the role of disorder is not clear. Although 
in the liquid phase the system is in a ballistic regime, 
when the conductance is quantized the disorder could 
take over in the localized phase and affect the charge 
distribution in the wire. AlAs wires, where the disorder 
is stronger, revealed conductance resonances explained in 
terms of Coulomb blockade (CB) physics^ CB behavior 
has also been found in the localized phase of GaAs wires.— 

Even if there are features that still need explanation, 
our calculations show that the electronic correlation plays 
a very important role at the experimental conditions, 
as the 2kF-to-Akp correlations transition takes place ex- 
actly in the proximity of the critical density for localiza- 
tion found in the experiment. In addition to this result, 
which is the main outcome of the paper, we also deter- 
mined the charge and spin velocities by means of the 
QMC method explained in Sec. [IT] and the effective J 
coupling via the WKB approach. We computed those 
quantities close to the transition for the homogeneous 
wire with r s — 1.25 (40/im -1 ). The charge velocity turns 
out to be v p = 2.33^f- The corresponding LL param- 
eter g = vf/v p = 0.43 is in agreement with previous 
estimates^ and comparable to the experimental g w 0.5, 
measured at the density of 40/im _1 F 7 - At the same den- 
sity the experimental value for vp jv a is in the range of 
1.1 — 1.6, while we found vpjv a — 1.24 at 40/j,to _1 . 

In Fig. [T2] we plot the full dependence of the spin ve- 
locities on the density computed with the perturbative 
generalized RPA (GRPA) ,22, WKB and the exact QMC 
methods. Although the GRPA is poor near the local- 
ization transition, it agrees with the QMC at high den- 
sity. As noted above the experimentally measured spin 
velocities are also in rough agreement with the QMC es- 
timate in a range of densities around n = 40/zm -1 . To 
show the importance of the microscopic details of the 
interaction in reproducing the measured values we also 
display in Fig. the GRPA prediction based on a dif- 



ferent model potential which assumes a screening due to 
a metallic gate at R = 50. 7 This latter model gives vir- 
tually unrenormalized spin velocities (v a w Vp) up to 
n = 40/jto~ 1 in contrast with the strong suppression of 
the values found in the experiment. 
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FIG. 12: (Color online) Inverse spin velocity of the in- 
finite wire. The red circles indicate estimates from the 
WKB approximation, whereas the black triangles are de- 
termined using the QMC method described in Section [TTJ 
The two lines are estimates due to the perturbative general- 
ized random phase approximation (GRPA)—: vf /v„ rpa — 
1 / \J\ — V(2kF)/(irvF)- The green line for the gated wire 
(with R = 50) uses the potential described by Auslaender et 
aZ.— whereas the dotted blue line uses the potential (Eq. I12p 
screened by the lower wire. 

Last but not least, our WKB estimate of J turns 
out to be of the order of the experimental temperature 
(T = 0.25-ftT) around n — 10/zm -1 . This means that at 
least the first few Coulomb blockade peaks in the exper- 
iment should be in a spin incoherent regime, where the 
LL description by Fiete et al. applies, although in the 
vicinity of the transition the spin degrees of freedom are 
not dominated by thermal broadening. 

V. CONCLUSIONS 

We have presented extensive quantum Monte Carlo 
calculations to study the properties of electrons con- 
strained to one dimension with a harmonic confinement 
and interacting via several different potentials. 

For unscreened interactions with a long range 1 jx tail 
there are three different regimes. At high density the 
electrons behave as a correlated liquid, transitioning to 
a quasi Wigner crystal as the density decreases, where 
strong 4kp correlations follow the LL predictions, 23 We 
accurately determined the crossover density for various 
thicknesses and found that the crossover is pushed to 
higher densities for thinner wires. Finally at very low 
densities the charge degrees of freedom are described by 
spinless fermions and the spins decouple with exponen- 
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tially small exchange interactions. We approached this 
limit by using the WKB approximation. 

When screening is introduced, the interactions are not 
long-range, and the quasi Wigner crystal order is de- 
stroyed. However, Ak p correlations are still present even 
in the case of screened interactions. The spinless fermion 
regime acquires a new behavior when the wire is very 
thin and the screening makes the potential short-range. 
In this case the particles act as though they were non- 
interacting and spinless in analogy to physics previously 
studied for bosons with an infinite contact repulsion. 26 

We applied our numerical approach to analyze a model 
chosen to realistically describe the double wire system 
studied in the experiments of Steinberg et al.fi- where 
a localization transition is observed. Our model assumes 
screening due to a second wire described within linear re- 
sponse theory, and includes the finite length of the wire 
via the external potential derived in Ref. |{J. We show 
that a crossover from a liquid to a state with Akp cor- 
relations occurs around the localization density found in 
the experiment. Additionally, our exact Monte Carlo cal- 
culations yield charge and spin velocities for this model 
in agreement with those observed in the experiment close 



to the transition. We stress that the observables such as 
the transition density and the spin velocity are partic- 
ularly sensitive to the microscopic details of the model 
interaction. To reproduce all features of the experiment 
it may be necessary to include further refinements such 
as a more accurate modulation of the external potential 
due to the gates, the effects of higher subbands in the 
transverse direction and the full treatment of interwire 
electronic correlation by explicitly including the electrons 
in the other wire. However, the simple model considered 
here shows that the exact treatment of electronic correla- 
tion is essential to quantitatively describe the localization 
transition seen in experiments. 
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